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Abstract. I discuss Haldane's concept of generalised exclusion statistics (Phys. 
Rev. Lett. 67, 937, 1991) and I show that it leads to inconsistencies in the 
calculation of the particle distribution that maximizes the partition function. These 
inconsistencies appear when mutual exclusion statistics is manifested between different 
subspecies of particles in the system. In order to eliminate these inconsistencies, I 
introduce new mutual exclusion statistics parameters, which are proportional to the 
dimension of the Hilbert sub-space on which they act. These new definitions lead 
to properly defined particle distributions and thermodynamic properties. In another 
paper (arXiv:0710.0728) I show that fractional exclusion statistics manifested in general 
systems with interaction have these, physically consistent, statistics parameters. 
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1. Introduction 

Haldane's concept of fractional exclusion statistics (FES) pQ have been applied to 
the study of many types of physical systems, revealing interesting properties. For 
example it has been applied to strongly interacting systems, such as Tomonaga- 
Luttinger model IH [5], Colagero-Sutherland model 0, El [9j [ID], fractional 
quantum Hall effect [HI [121 or to interacting particles in one or two-dimensional 
systems, described in the mean field approximation [HI [El EH EE]- The statistical 
properties of FES systems have been calculated mainly by Isakov [18] and Wu [19J, 
while Iguchi extended the Fermi liquid model to the model of a FES liquid [201 ET4 ; 
the microscopic reason for the manifestation of FES have also been discussed by several 

authors nn ens nni 121 cna 13 ng. 

Although the concept received so much attention and has been applied in general 
to many types of systems, I will show here that when mutual exclusion statistics 
is manifested between different subspecies of particles in the system, FES leads to 
thermodynamic inconsistencies. I will also show that these inconsistencies can be 
corrected by a redefinition of the exclusion statistics parameters. 

In a related paper I showed that fractional exclusion statistics appears in general 
systems of interacting particles and the statistics parameters indeed obey the rules 
conjectured here [23] . 

2. Thermodynamic inconsistencies in FES 

In this section I will prove using two model systems, that in FES systems the equilibrium 
particle populations are ambiguously defined, if mutual statistics parameters are not 
zero. For this, I will recalculate the partition function and the most probable particle 
distribution in a FES system, following the procedure used by Wu in Ref. [19]. 

Haldane defined the fractional exclusion statistics as acting on Hilbert spaces of 
finite dimensions p]. If we have only one such a space, in which we put N ideal bosons 
or fermions, then the number of microscopic configurations we have in the system is 
W b = (G + N - 1)!/[JV!(G? - 1)!] (for bosons) or W f = G\/[N\(G - N)\] (for fermions). 
Fractional exclusion statistics of parameter a is an interpolation between these two cases 
and the number of configurations is W = [G + (N - 1)(1 -a)]\/{N\[G-aN - (1 -«)]}- 
we say that the addition of 5N particles in the system reduces the number of available 
states in the system by aSN [TJ [19]. 

Now let us generalize the problem to the case when we have more than one Hilbert 
space. The spaces are Hq, 7ii, . . . , of dimensions G , G\, . . . , and which contain N , Ni, 
. . . , particles. In this case we have the FES parameters with i, j = 0, 1, . . .. Mutual 
exclusion statistics is manifested between the spaces Hi and Hj {i ^ j) if ^ 0-we say 
that the addition of SNj particles in the space Tij changes the number of available states 
in the space Tii by —aij5Nj. With these notations, the total number of configurations 
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is 



w = u 



[Gi + Ni-l-^jotijiNj-Si. 



NA 



G, 



(1) 



1 - EjCtijiNj - Sij 

Having the number of microscopic configurations ((Tj), if we asign the energy and the 
chemical potential /ij to the states in the space i, we can calculate the grandcanonical 
partition function, Z 



(2) 



Z= J2W({N t })exp 
W} 

and the total energy of the system in the given configuration, E = J2i -ZVjej-we use the 
notation (3 = l/ksT, where T is the temperature of the system. 

The most probable configuration, {Ni}, is obtained by maximizing Z with respect 
to the set {Ni}. If we introduce the notations Hi = Ni/Gi and = a^Gj/Gi, and 
assume that for each i both, Gi and Ni are sufficiently large, so that we can use the 
Stirling approximation [lnGj! ~ Gi \n(Gi/e) and lniVj! ~ Ni ln(A^/e)] the maximization 
procedure gives us the system of equations, 

1 — Sfc ftjk n k 



me 



1 + XX<^ _ Pik)nk 



n 



1 + Efe(^ife - Pjk)n k 



(3) 



The system (J3j) is solved easier if we denote Wi 
(IHl) becomes 



[1 + Wi 




■J2k Pik n k/ n i- Using this notations 



(4) 



and niS can be calculated from the new system, 
^(SijWj + fcj)nj = 1. 



(5) 



If the spaces i are fixed, Eq. (jSJ) can be solved and the populations rii can be 
determined. But in general, the spaces Hi can be changed by dividing or combining 
them and, if the system has a proper thermodynamic behavior, such changes should not 
affect the thermodynamic results. Unfortunately this is not the case and I will show now 
by simple examples that the particle distribution that maximizes the partition function 
depends on our choice of subspaces TCi. 

Let us assume that we have only two finite dimensional Hilbert spaces, 7io and Tlx, 
of dimensions Go and G\, and the FES parameters aoo, atn, ®oi, and a\Q. Then Eqs. 
(j3J) and (jl]) reduce to systems of two equations and two unknowns: 

- g/3(eo-Mo) 

(6) 



1+WQ 



"01 



«11 
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and 



(w + a 00 )n + a i^^i 



«io§7% + (wi + a n )ni 



(7) 



1 



respectively. Now let us further assume that the space Hi was obtained as a union of the 
two smaller, disjoint subspaces, Hi a and Hu, of dimensions G\ a and Gib, respectively. 
Then, instead of two Hilbert spaces and four exclusion statistics parameters that we 
started with, we can also describe the system as consisting of three subspaces and 
nine statistics parameters. From these nine, «oo remains unchanged, whereas otoi a 
and cioib should be identical to aoi- A natural choice for the other six parameters 
is a lala = a lblb = an, a lb0 = a la0 = a w , and a lalb = a uia = 0. Obviously, we also 
have 6i a = exb = ei and ji\ a = = pi- In the new configuration, the systems (jlj) and 
(15]) become 











) 







D /3(<E -Mo) 



, an 



3 /9(ei-A»i) 



(8) 



and 



(w' + a 00 )n + a 01 ^n la + a 01 ^n lb = 1 



(9) 



, "locft^o + (w u + a u )n lb = 1 

The systems ([6]) and (jTJ) should admit the same physical solution as the systems (jHJ) 
and ([9]), which is n\ a = n\ b = n\. If we also take, for the simplicity of the calculations, 
G\ a = Gib = Gi/2, then from the last two equations of Qj we get wi a = wi b and 

a w n G /Gi + (w la - wi)ni = 0. (10) 

Using wi a = Wu = w' 1 into (jSJ), we obtain the system 

, «oo / ,„' \ 2cno 



o f3(e -fj, ) 



(11) 



Since the systems flS} and (TTTT) are not identical, Eqs. (TTTT) and (ITU|) (eventually combined 
also with [6]) give three equations for the two unknowns, w' and u^. The system of 
equations is overdetermined and does not have solutions, unless the statistics parameters 
have some special values. Such special values are e.g. «io = «oi — («oi — would be 
necessary if we would split the space Ho instead of Hi). 
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2.1. Large systems 

I will show now that the thermodynamic inconsistency of FES happends also in systems 
of infinite, quasi-continuous spectra. For this let us assign to the single-particle states 
the sets of quantum numbers, k-the set k does not necessary consist of wave-vector 
components, as it usually does, but it may contain any quantum numbers. If the system 
is large enough, we say that each of the subspaces Hi contains the set of single particle 
states {k}j, which "fill-up" a volume Vi in the space of the quantum numbers and we 
define the density of states, a^, so that Gi = f v . dk and iV, = J v _ n k o"k dk. Then 
the parameters depend on the quantum numbers {k}j and {k}j, so we shall change 
the subscripts and write in general, e.g. akk'- In these new notations, the system (Tj0) 
becomes 



for any division of the total Hilbert space of the system into elementary volumes Vi, 
as long as Gi and Ni (for any i) are large enough, so that we can apply the Stirling 
formula and the maximization procedure presented above. Moreover, if we can choose 
all the volumes Vi so that they are small enough to use the approximation Gi = a^Vi 
and c*ij = c*kk' for any k in Vi and k' in Vj, Eq. (jSJ) becomes 



where we also changed the notation Vi into Vk and Vj into V\j. Now we observe directly 
that while Eqs. ffl2|) contain only intensive parameters (i.e. which do not depend on the 
volumes Vk involved), the values of the equilibrium particle populations f|T3|) depend 
on our arbitrary choice of volumes, so the thermodynamic quantities, which are the 
populations rik, are not well defined. 

3. Correction of the parameters 

The reason for which the thermodynamics of FES systems is ambiguous is obvious from 
the beginning. If we look for example at Eq. ([3]) and we imagine that we reduce the 
dimension of one of the spaces, say of space i, by half, than since all the parameters 
Pij are proportional to Gj/Gi, this could cause a significant resistribution of all the 
occupation numbers (as we proved in concrete cases above). Since the population 
numbers should be intensive parameters, it means that the system could not have 
thermodynamic limit. 

Fortunately, as it was easy to observe the anomaly of FES, it is as easy to correct it. 
I shall prove below that if the "trouble-making" mutual exclusion statistics parameters 
are proportional to the space on which they act, i.e. oc Gi, then the system admits 
a thermodynamic limit. In Ref. [23] I analized the exclusion statistics that appears in 
general interacting systems and I showed that in the limit of a quasi-continuous density 
of states the mutual exclusion statistics parameters indeed satisfy this property. To 




(12) 




(13) 
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differentiate these new parameters from the ones used above, I shall denote them by 
<Skk'- Therefore, the qualitatively new feature of the parameters cikk' is that for k 7^ k', 
«kk' is not a simple number, but it is proportional to the dimension of the subspace V k : 
«kk' = GV k «kk' = CkVk^kk'- If k = k', the statistics parameters are similar to the ones 
studied so far and we write <Skk = ctkk- 

Let us now try to find the equilibrium configuration with these new statistics 
parameters and see if this has a thermodynamic limit. To do this, we first have to 
write down the new number of configurations (CQ): 



w = \{ 



[G Vk + (1 - a kk )(AV k - 1) - G Vk Ek^ k akk'iV v 



NA 



Gy. - 1 



(14) 



"kk(AV k - 1) - G Vk Ek'^ k ttkk'A^y k , 

But we notice already that for large enough volumes V k , both, Gy k and Ny k are large, 
therefore the summations Gy k J2 k >^ k a kk' Ny k , , containing terms bilinear in these large 
quantities, must be much larger than Gy k and Ny k . In conclusion both square brackets 
in ([141) should be negative. This is not admissible. [We note in passing that also in 
Haldane-Wu formalism (TjQ) [TJ [19] negative terms may appear in the calculation of the 
number of configurations for unappropriate choice of subspaces, even though in the 
summations there are only terms linear in iVj.] 

To solve the problem we have to regard the exclusion statistics as acting on the 
quasiparticle added as a perturbation to the equilibrium quasiparticle distribution. 
Therefore we assume that on top of the equilibrium distribution of particles, say N Vk , 
we add a small perturbation, 5Ny k , which changes the number of configurations into 



n 



[Gv k + N Vk + 6N Vk - 1 



Ek^ k «kk'£AT v 



(N Vk + 5N Vk )\(G Vk -l-E 



k ^ k a kk >5N v ] 



(15) 



[G Vk + N Vk + (1 - a kk )5N Vk -l-G Vk Ek^ k cWiV v 



(N Vk + SN Vk )\ 



G Vk 



a kk 5N } 



Vk 



Gv k Sk'^ k cx kk >5N v 



Notice that now Gy k is the Bose dimension of the Hilbert space, i.e. the number of 
available states, and not its real dimension. So, unless we have an ideal Bose system, 
the ratio Ny k /Gv k , denoted by h k , is actually the ratio of the number of particles 
to the number of available states (or holes). In terms of G k of equation (pQ) or (CO 



Gy k = G k — aN k . The equilibrium distribution is then obtained by imposing that the 
partition function is stationary with respect to the perturbations SNy k . If we introduce 
the notations, 5h k = 5Ny k /Gv k , using the Stirling approximation we write the logarithm 
of W as 



1 + ™k + (1 - a kk )5n k - MKk'Gv^Sni 



x In 



k'(^k) 



1 + n k + (1 - a kk )5h k - a kk>G Vk ,Sn k/ 

k'(^k) 



- (n k + 5n k ) ln(n k + 5n k ) - ( 1 - a kk <fa k - 1 - Y a *-k'G Vk ,Sn k > 

k'(^k) 
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x In 1 - a kk <5ri k - 1 - ^ «kk'<5k'^k' \ (16) 
V k'(^k) / J 

Adding (3 J2k GV k n k (/i k — £k) to (|T6|) . we obtain the logarithm of the partition function, 
which, if maximized, gives 

d 







In W + P G v k , n k' (pv ~ e k ' 



= G Vk I (3(nk ~ £k) - hin k + (1 - ct kk ) ln(l + n k ) - XI "k'kGV k , hi(l + n k /) > (17) 

I k'(^k) J 

We observe now that if ct kk / = 0, for any k and k', we obtain the Bose distribution, 
nk = nk — {exp[/3(e k — /i k ) — whereas if a kk / = 0, for any k ^ k', and a kk = 1 

for any k, we obtain n k = exp[/3(// k — e k )], which is the ratio of the number of particles 
to the number of holes in a Fermi system, as expected. 

Now we can finally extract from (fTTj) a self-consistent equation for n k : 



P(Vk-ek) + \n [ -—-^ = ]T Gy k/ ln[l+n k ,]a k , k (18) 

nk k'(^k) 

In the thermodynamic limit, when Gy k <C J2k'(^k) Gv k n Eq. (Tl8|) gives a good, 
unambiguous particle distribution in the system. The summation to the right-hand 
side of Eq. (fl8l) can be readily and consistently transformed into an intergral when the 
dimension of the system increases and we get the self-consistent integral equation for 
the particle distribution, 

P(»k - e k ) + In l - ^ = / a v ln[l + n k ,]a k , k dk'. (19) 

n k J 

In the end let us remark an identity. If all the mutual statistics parameters were 
zero, the model introduced here should be identical to the previous FES model (section 
[2]) so so let us compare the results. First, if a kk / = for all k ^ k', equation ffl8l) 
reduces to 

(l + ^k) 1 — 



e^" k \ (20) 
n k 

whereas equation ( fl2|) becomes 

„,"kk 

^ = e^~^ (21) 

(l+W k ) Q kk-l V ' 

Now one can check easily that equation ( 1201) is identical to equation ( 12T1) if n k = . 
But on the other hand, by definition, w^ 1 = N k /(Gk — a kk A^" k ) = N^/G^ = n k . 
Therefore the results are indeed identical and we observe with this ocasion that in 
FES systems with no mutual statistics, the quantity w k was the ration between the 
Bose dimension of the subspace Vk and the number of particles in this subspace. 
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4. Conclusions 

I showed that fractional exclusion statistics leads to ambiguous thermodynamic results 
if the mutual exclusion statistics parameters are not zero (a^ 7^ 0, for i 7^ j). To correct 
this ambigiuty and to obtain consistent thermodynamic results I introduced new mutual 
exclusion statistics parameters, which are proportional to the dimension of the subspace 
on which they act. 

In a related publication [23] I proved that a gas of interacting particles can be 
described as a fractional exclusion statistics gas, with the exclusion statistics parameters 
having the properties obtained here. 

Acknowledgments 

I thank Dr. A. Parvan for motivating discussions. This work was done at the Bogoliubov 
Laboratory of Theoretical Physics, JINR Dubna, Russia and I thank the stuff of the 
Laboratory, especially Dr. S. N. Ershov and Dr. A. Parvan, for hospitality. The work 
was partially supported by the NATO grant, EAP.RIG 982080. 

References 

[1] F. D. M. Haldane. Phys. Rev. Lett., 67:937, 1991. 

[2] S. Tomonaga. Progr. Theor. Phys., 5:544, 1950. 

[3] J. M. Luttinger. J. Math. Phys., 4:1154, 1963. 

[4] D. C. Mattis and E Lieb. J. Math. Phys., 6:304, 1965. 

[5] J. M. P. Carmelo, P. Horsch, A. A. Ovchinnikov, D. K. Campbell, A. H. Castro Neto, and N. M. R. 
Peres. Phys. Rev. Lett., 81:489, 1998. 

[6] F. Colagero. J. Math. Phys. (N.Y.), 10:2191, 1969. 

[7] B. Sutherland. J. Math. Phys. (N.Y.), 12:247, 1971. 

[8] B. Sutherland. Phys. Rev. A, 4:2019, 1971. 

[9] B. Sutherland. Phys. Rev. A, 5:1372, 1972. 

[10] M. V. N. Murthy and R. Shankar. Phys. Rev. B, 60:6517, 1999. 

[11] A. D. de Veigy and S. Ouvry. Phys. Rev. Lett, 72:600, 1994. 

[12] T.H. Hansson, J.M. Leinaas, and S. Viefers. Nucl. Phys. B, 470:291, 1996. 

[13] S.B. Isakov and S. Viefers. Int. J. Mod. Phys. A, 12:1895, 1997. 

[14] M. V. N. Murthy and R. Shankar. Phys. Rev. Lett., 73:3331, 1994. 

[15] D. Sen and R. K. Bhaduri. Phys. Rev. Lett., 74:3912, 1995. 

[16] R. K. Bhaduri, S. M. Reimann, S. Viefers, A. G. Choudhury, and M. K. Srivastava. J. Phys. B, 
33:3895-3903, 2000. 

[17] T. H. Hansson, J. M. Leinaas, and S. Viefers. Phys. Rev. Lett., 86:2930-2933, 2001. 

[18] S. B. Isakov. Phys. Rev. Lett., 73(16):2150, 1994. 

[19] Yong-Shi Wu. Phys. Rev. Lett, 73:922, 1994. 

[20] K. Iguchi. Phys. Rev. Lett., 80:1698, 1998. 

[21] K. Iguchi. Phys. Rev. B, 61:12757, 2000. 

[22] K. Iguchi and B. Sutherland. Phys. Rev. Lett, 85:2781, 2000. 

[23] D. V. Anghel. larXiv:0710.0728l 2007. 



